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ABSTRACT 

The Balnicr line profiles of nonradiative supernova remnant shocks provide the means to 
measure the post-shock proton velocity distribution. While most analyses assume a Maxwellian 
velocity distribution, this is unlikely to be correct. In particular, neutral atoms that pass through 
the shock and become ionized downstream form a nonthermal distribution similar to that of 
pickup ions in the solar wind. We predict the Ha line profiles from the combination of pickup 
protons and the ordinary shocked protons, and we consider the extent to which this distribution 
could affect the shock parameters derived from Ha profiles. The Maxwellian assumption could 
lead to an underestimate of shock speed by up to about 15%. The isotropization of the pickup ion 
population generates wave energy, and we find that for the most favorable parameters this energy 
could significantly heat the thermal particles. Sufficiently accurate profiles could constrain the 
strength and direction of the magnetic field in the shocked plasma, and we discuss the distortions 
from a Gaussian profile to be expected in Tycho's supernova remnant. 

Subject headings: supernova remnants-sliock waves-ISM: lines and bands-line: profiles-turbulence 



1. Introduction 

Fast interstellar shock waves that encounter 
partially neutral gas are observable as filaments of 
pure Balmer line emission if they are young com- 
pared to their radiative cooling times (Chevalier 
& Raymond 1978; Raymond 1991; Ghavamian et 
al. 2001). The Balmer lines are produced in the 
thin layer just behind the shock where hydrogen 
atoms are excited and ionized, and this layer is 
thin enough that Coulomb collisions cannot bring 
different particle species into thermal equilibrium. 
Hence the Balmer lines can be used to probe the 
physical processes in coUisionless shocks. 

The Balmer lines have two component line pro- 



files. The broad component arises from neutral 
H atoms created by charge transfer with post- 
shock protons, and its velocity width is compa- 
rable to the downstream proton thermal velocity. 
The narrow component comes from neutrals that 
have passed through the shock, but that have 
not been ionized by charge transfer. Therefore, 
its velocity width corresponds to the tempera- 
ture of the pre-shock gas. The intensity ratio of 
the broad and narrow components depends on 
electron and proton temperatures, Tg and Tp, so 
that it can serve as a diagnostic for T^/Tp im- 
mediately behind the shock. This is an impor- 
tant quantity for interpreting X-ray spectra of 
SNRs and for understanding coUisionless shocks. 
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In a few cases it has been possible to measure 
the widths of UV hnes of other elements, and 
theref ore the kinetic tempe r atures, Tj, of other 
ions (Raymond et al. 1995 : Laming et al. 19961 : 



Raymond et al. 2003t iKorreck et al 20071) . The 
overall result is that the plasma behind rela- 
tiyely slow shocks 300 km s~^ ) is close to 
thermal equilibrium, while shocks faster than 
about 1000 km s~^ are far from e quilibrium, with 
Te/r„ < 0.1 and T,/T„^ rm/mp (|Rakowski 20051: 
Ghayamian et al. 20071 ). Other important appli- 
cations of Balmer line diagnostics for coUisionless 
shocks are estimates of shock speed, which can be 
combined with proper motions to find SNR dis- 
tances (Winkler, Gupta & Long 2003), and infer- 
ences of cosmic ray diffusion c oefficients from the 
properties of shock precursors (ISniith et al. 1994 : 
Hester et al. 1994HLee et al. 20071 ). 

All of the current models used to interpret the 
Balmer line profiles assume that the post-shock 
proton yelocity distribution is Maxwellian (Cheya- 
lier et al. 1980; Lim & Raga 1996; Laming et al. 
1996), though there is no solid justification for that 
assumption. Coulomb collisions are not able to 
bring the protons to a Maxwellian rapidly enough, 
and it is not clear what sort of distribution would 
be produced by plasma turbulence. Heng & Mc- 
Cray (2007) have recently drawn attention to the 
importance of sequential charge transfer events in 
determining the profile of the broad component at 
shock speeds above about 2000 km s^^, where the 
charge transfer cross section changes rapidly with 
energy, and this affects some of the diagnostics. 
Heng et al. (2007) have extended the model effort 
to a fuller treatment of the hydrodynamics than is 
usually employed, but they keep the assumption 
that the proton distribution is Maxwellian. 

Ion yelocity distributions directly measured in 
the solar wind are essentially never Maxwellian in 
the vicinity of shocks. Ion distributions typically 
have power-law tails or strong anisotropics, with 
beamed components upstream and highly perpen- 
dicular enhancements downstream (Schopke et al. 
1983; 1990; Goshng and Robson 1985; Thomsen 
1985; Kucharek et al. 2004). The Balmer hue 
profiles of non-radiative shocks provide a unique 
opportunity to search for non-Maxwellian yelocity 
distributions in astrophysical plasmas. 

In this paper we will keep the assumption 
that protons that pass through the shock have a 



Maxwellian distribution at the temperature given 
by the Rankine-Hugoniot jump conditions (Draine 
and McKee 1993), but we will add the manifestly 
non-Maxwellian distribution of protons that pass 
through the shock as neutrals and become ion- 
ized. We will consider the potentially observable 
effects on the Balmer line profiles including line 
widths and centroid shifts and how they might 
affect shock parameters derived from Ha profiles. 
We will also discuss the implications of magnetic 
field strength and direction and of plasma tur- 
bulence on the profiles and the possibility that 
observed profiles could constrain the field param- 
eters. We will briefly consider the implications 
of non-Maxwellian distributions for the heating of 
electrons and minor ions. 

2. Pickup Ions 

Neutral particles that are ionized in the post- 
shock flow are very much like the pickup ions 

(PUI) measured by spacecraf t in the solar wind 

iMoebius et al. 19851 : . Gloeckler et al. 19931 : llsenberg 1995f) . 
When neutral atoms slowly flow into the inter- 
planetary medium, they can be ionized by pho- 
tons from the Sun, by charge transfer with solar 
wind ions, or by collisions with electrons. At 
that point, the new ions are streaming with re- 
spect to the solar wind plasma at the solar wind 
speed, Vsw: which is much larger than the local 
Alfven speed, Va- These ions are immediately 
swept up by the magnetic field in the solar wind. 
Their velocity component perpendicular to the lo- 
cal magnetic field becomes a gyrovelocity around 
the field, which, in combination with the instan- 
taneous parallel component, initially forms a mo- 
noenergetic ring-beam in velocity space. This 
ring-beam is unstable, and the particles rapidly 
scatter toward isotropy by interacting with ambi- 
ent or self-generated waves, resulti ng in a velocity- 
space shell (Sagdeev et al. 1986; iLee fc Ip 1987 : 
Isenberg 2005i : feogdan et al. 199ll ). 

In the solar wind, pickup protons are distin- 
guished by their unusual velocity distributions, 
but heavier pickup ions can also be recognized 
by their single ionization states, such as He"*" or 
O"*", which stand out among solar wind ions that 
are much more highly ionized. The pickup ions 
add significant mass and momentum to the so- 
lar wind in the outer heliosphere, and the waves 
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they generate play an import ant role in heating 



the solar wind beyond 5 AU ([Smith et al. 2001 
llsenberg et al. 20031 : llsenberg 200^ . The com 



position and charge state of these pickup ions 
indicate that they form the source particles for 
the observed anomalous component of cosmic rays 
(Garcia-Munoz et al. 1973, 1975; Fisk et al. 1974; 
Cummings & Stone, 2007). These particles must 
therefore be preferentially accelerated to several 
tens of MeV/nucleon at the solar wind termination 
shock or in the heliosheath beyond. However, re- 
cent in situ observations during and after the Voy- 
ager encounters with the termination shock (Stone 
et al. 2005; Stone 2007) have shown that the en- 
ergization process is still not well understood. 

3. Consequences for SNR shocks 

Consider a planar shock in which the down- 
stream magnetic field makes an angle 9 with the 
shock normal. Since the field component perpen- 
dicular to the flow is compressed by the shock, 
is typically 60° to 85°, though of course pure 
parallel and pure perpendicular shocks maintain 
their field directions. For a strong shock with a 
compression ratio of 4, a neutral passing through 
the shock moves at jVs relative to the post-shock 
protons. Thus when it becomes a pickup ion it 
acquires a gyro velocity 



Vi_ = -Vs sin{9) (1) 
and a velocity along the field direction of 



Vi\ = -Vs cos{e) (2) 

relative to the post-shock plasma. These monoen- 
ergetic particles form an unstable ring distribution 
in velocity space. They can emit plasma waves 
and interact with these waves to scatter into a 
more isotropic distribution. Generally, the dom- 
inant isotropization process is pitch-angle scat- 
tering through the cyclotron resonant interaction 
with parallel-propagating ion-cyclotron and fast- 
mode waves (Wu & Davidson 1972; Winske et al. 
1985; Lee & Ip 1987; see also Zank 1999; Szego et 
al. 2000). 

The ring-beam distribution may also be subject 
to other plasma instabilities, depending on the rel- 
ative density and downstream conditions. In prin- 



ciple, a downstream magnetic field nearly paral- 
lel to the flow can result in bump-on-tail (Gary 
1978) or firehose-like instabilities (Winske et al. 
1985; Sagdeev et al. 1986). The saturation of the 
Landau bump-on-tail instability leaves a highly 
anisotropic beam which still scatters in pitch angle 
through the cyclotron resonance. The firehose in- 
stability could disrupt the beam, but requires both 
a high density of pickup ions relative to the back- 
ground plasma and an ionization time-scale much 
shorter than the time-scale for cyclotron resonant 
pitch-angle scattering. If the downstream mag- 
netic field is nearly perpendicular to the flow, a 
mirror-mode instability can be excited (Winske & 
Quest 1988; McKean et al. 1995), but this insta- 
bility saturates at a much lower level than the reso- 
nant ion-cyclotron instability (Yoon 1992), and so 
may be neglected. In this paper, we will take the 
ring-beam of newly-ionized protons to quickly sta- 
bilize through cyclotron-resonant pitch-angle scat- 
tering. In particular, we will assume the rapid 
formation of a bispherical distribution. 

3.1. Bispherical Distribution 

Under most conditions, a given energetic 
proton is cyclotron-resonant with two parallel- 
propagating electromagnetic modes. If the pro- 
ton parallel speed is much faster than the Alfven 
speed, Va, both of these waves will be Alfven 
waves - one propagating along the field in the 
same direction as the proton and the other in 
the opposite direction. Resonant scattering away 
from the ring-beam will result in the amplifica- 
tion of one of these modes and the damping of 
the other. Which mode is unstable depends on 
the position of the ring-beam in velocity space, as 
determined by the angle of the local magnetic field 
to the plasma flow direction. The resonant inter- 
action with either wave yields a diffusion which 
conserves the proton energy in the frame of the 
wave phase speed, scattering the particles along 
a sphere in velocity space centered on one of the 
points W|| — ±Va, as shown in Figure 1. A useful 
analytical result is obtained in the case where the 
damped mode can be neglected and the scattering 
at each point in velocity space is only due to in- 
teractions with the unstable mode. In this case, a 
steady ring-beam will be scattered to a bispherical 
distribution: a uniformly populated shell formed 
by the two spherical caps which meet at the posi- 
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tion of the original ring-beam (Galeev & Sagdeev 
1988; Williams & Zank 1994). 

Many of the basic properties of this distribution 
may be obtained geometrically. If the ring-beam 
of the newly ionized protons is located at {V\\, V±) 
as given by equations (1) - (2), the radii of the two 
spherical caps are Vj_ — VI + {V\\±Va)^- The area 
of each cap in velocity space is a± = 2ttv±{v± =f 
V|| — Va). Since the particles are distributed uni- 
formly over these areas, the net streaming speed 
of the bispherical shell is 



Vbuik = — [Vk(a+ - a-) + ttV^{v^ - v+)] (3) 
ax 

where the total shell area = a+ + a_. Clearly, 
the case of flow perpendicular to the magnetic 
field, V|| = 0, gives v_^. = w_ and v^uik — 0. Simi- 
larly, for parallel flow faster than the Alfven speed, 
the distribution reduces to a single sphere of ra- 
dius — Va, and the bulk speed is slowed to 
Vbuik = Va- In general, the streaming speed of 
the bispherical distribution is bounded by ±Va- 
Figure 2 shows this streaming speed as a function 
of magnetic field angle 9 for several values of the 
downstream field strength, taking a shock speed 
of 2000 km and an upstream proton density 
of 1 cm~'^. 

These simple properties may be modified for 
realistic conditions. For instance, dispersion of 
the resonant waves will systematically shift their 
phase speed, and so distort the shape of the shell 
away from a sphere (Isenberg & Lee 1996). This 
distortion can be significant if the speed differ- 
ence between the neutrals and the downstream 
plasma is comparable to Va- In addition, an effi- 
cient turbulent cascade could maintain the stable 
wave mode intensity despite the damping by the 
pickup protons. In this case, the multiple wave- 
particle interactions with both stable and unsta- 
ble waves can yield a much different distribution, 
and even result in particle acceleration through 
the second-order Fermi mechanism (Isenberg et al. 
2003; Isenberg 2005). However, the bispherical ex- 
pressions provide a reasonable first approximation 
to the pickup proton distribution expected down- 
stream of a strong supernova shock. In this initial 
study, we will retain the bispherical assumptions, 
and address these simplifications in the discussion 
section. 



3.2. Total Proton Distribution 

At any point in the downstream plasma the ve- 
locity distribution is the sum of the distributions 
of the shocked protons and the protons formed by 
ionization or charge transfer in the downstream 
gas. If the preshock neutral fraction is small, 
the distribution is dominated by the shocked pro- 
tons, and it will be difficult to detect the ef- 
fects of the pickup protons. These effects will 
be much easier to see in the shocks of Tycho's 
SNR, where the neutral fraction is around 0.85 



(jGhavamian et al. 20001) than in SN100 6. where 
it is around 0.1 (|Ghavamian et al. 2002h . Figure 
3 shows a simple model of a shock propagating 
at 2000 km into a medium with uh = rip = 
0.5cm~^, roughly similar to the values expected 
for Tycho's SNR. The proton density just behind 
the shock is the density of thermal protons, so 
the increase downstream represents the addition of 
pickup protons. The neutrals immediately behind 
the shock make up the slow or narrow component. 
Their density drops as charge transfer converts 
them to pickup protons and relaces them with fast 
or broad component neutrals. Eventually, colli- 
sional ionization removes all neutrals, leaving a 
fully ionized plasma far downstream. The rate 
coefficients for charge transfer and ionization by 
electrons and protons were adopted from Lam- 
ing et al. (1996). Note that this plot assumes 
that the pickup protons move with the same bulk 
speed as the background plasma. This will be 
strictly true only for a perpendicular shock, since 
the scattered shell of pickup ions will generally re- 
tain some streaming motion with respect to the 
thermal plasma if the field has a component along 
the flow. 

Figure 4 shows the thermal proton and pickup 
ion distributions for one choice of the parame- 
ters. For the modest Alfven speeds expected be- 
hind SNR shocks, the PUI distribution is not far 
from spherical. Thus the total velocity distribu- 
tion shows a peak with a sharp cutoff plus high 
velocity wings from the thermal distribution. 

The broad components of the Balmer line pro- 
files will reflect the proton distributions, though 
they are weighted by the charge transfer cross sec- 
tion. Figure 5 shows the proton velocity distribu- 
tion profiles in the direction parallel to the shock 
front obtained by adding the background plasma 
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distribution to the PUI distribution. Wo have cho- 
sen this direction because strong limb brightening 
is required to make the Ha emission from a non- 
radiative shock bright enough that a high S /N pro- 
file can be obtained. The projection is obtained by 
multiplying the velocity along the magnetic field 
direction by sm6. If the observer is not in the plane 
containing the pre-shock and post-shock magnetic 
field, the centroid shift will be reduced. In this pa- 
per we do not compute Balmer line profiles, since 
they depend on specific shock parameters. Such 
calculations will be needed for the interpretation 
of observations, but for shocks below roughly 2000 
km s~^ the variation of charge transfer cross sec- 
tion with velocity is weak enough that the Ha 
profile should closely resemble the proton veloc- 
ity distribution (see Heng and McCray 2006). It 
should be kept in mind, however, that the broad 
component of Ha is emitted from a region of vary- 
ing pickup ion fraction, with values near zero near 
the shock and approaching the pre-shock neutral 
fraction far downstream. Roughly speaking, the 
Ha profile will correspond to a pickup ion fraction 
of about half the pre-shock neutral fraction. 

It is apparent from Figure 5 that the departure 
from a Maxwellian ought to be detectable with 
sufficiently high S/N data. The difficulty is that 
the narrow component, whose intensity is gener- 
ally dominant, obscures the center of the broad 
component profile. The usual procedure of fitting 
the sum of two Gaussians to the total profile pro- 
vides enough degrees of freedom to absorb modest 
departures from the assumed Gaussians, especially 
if the far wings of the profile and the background 
level are poorly defined. For very fast shocks, the 
dropoff in charge transfer cross section with ve- 
locity may suppress the high velocity tail in any 
case. 

As an estimate of the error that could be made 
by assuming a Maxwellian proton distribution and 
using the resulting broad line width to derive a 
shock speed, we fit the profiles in Figure 5 with 
single Gaussians and compared those widths to 
the widths of a pure Gaussian at the tempera- 
ture expected from shock speed. We find that 
the Gaussian widths estimated from the Figure 
5 distributions are as much as 14% narrower than 
those predicted for a pure thermal distribution of 
protons, so the shock speed could be underesti- 
mated by 14%. This is the extreme case, however, 



and underestimates about half that large would 
be typical. These underestimates would be partly 
countered if the pickup process provides additional 
heating to the plasma. 

3.3. Plasma Heating 

Another possible consequence of the pickup 

process is plasma heating by the waves generated 
in the isotropization of the initial ring-beam of 
newly ionized protons. The energy lost by the 
protons in scattering from the ring-beam to the fi- 
nal nearly isotropic shell is transferred to the res- 
onant waves. These waves in turn may heat the 
plasma, either directly or through a turbulent cas- 
cade to dissipative modes. In the simple bispher- 
ical picture of section 3.1, the energy available to 
the waves is given by 

Eyj = Ec, — Ebd+ — Ebd- (4) 

where Eo = mn{V]_ -\- V^)/2 is the energy in the 
initial ring-beam and the energy in the bispherical 
distribution is 

Ebd± = ':^^[\v\^VA±ViTvl)+{vi-VAr] 

(5) 

Figure 6 shows the ratio of the total bispherical 
energy, Ebd = Ebd++Ebd- to the initial energy 
for various combinations of the Alfvcn speed and 
the downstream magnetic field angle. The wave 
energy in (4) is essentially a maximum estimate, 
since the; bispherical distribution has a lower en- 
ergy than the distributions obtained by including 
dispersive effects or the replenishment of the sta- 
ble wave modes (Isenberg & Lee 1996; Isenberg 
2005). 

The form of the plasma heating which results 
from the pickup proton generated waves is an ac- 
tive area of research in the solar wind. A phe- 
nomenological model which assumes that these 
waves feed a turbulent cascade which dissipates 
by heating the thermal ions has been shown to 
reproduce the observed proton temperatures in 
the outer heliosphere reasonably well (Smith et 
al.2001, 2006; Isenberg et al. 2003; Isenberg 2005). 
This heating can be important when the upstream 
neutral fraction is large, and it may therefore affect 
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the estimates of the shock speed from the observed 
B aimer hne width. 

3.4. Electron Heating 

Electron heating is observed to be very ineffi- 
cient in fast shocks, so the observed electron tem- 
peratures could provide a strong constraint on the 
wave energy even if only a modest fraction of the 
wave energy is transferred to the electrons. Ob- 
servations of young SNRs show that T^/Ti is less 
than 0.1 in shocks faster than abont 1500 km/s 
(Rakowski 2005). Ghavamian et al. (2006) pro- 
pose that cosmic rays diffusing ahead of a fast 
shock produce lower hybrid waves which then heat 
the electrons to a temperature of about 0.3 keV, 
and this can reproduce the observed variation of 
Te/Tj with shock speed. 

Alternatively, if the pickup proton ring distri- 
bution generates lower hybrid waves, they could 
heat electrons. The lower hybrid heating is ineffi- 
cient unless the Alfven speed is large (Omelchenko 
et al. 1989; Cairns fc Zank 2002), but the Alfven 
speed downstream of SNR shocks is very poorly 
known. In the absence of information about Va, 
one cannot make quantitative predictions. In Ty- 
cho's SNR, which has a large neutral fraction in 
the pre-shock gas, the observed low electron tem- 
perature precludes efficient transfer of energy to 
the electrons if Va > O.lVg. 

3.5. Downstrecim Heating of Heavy Ions 

Other elements with ionization potentials at 
least as large as that of hydrogen will be partially 
neutral in the pre-shock gas. They will also be ion- 
ized and picked up except that they will be more 
likely to undergo electron or proton coUisional ion- 
ization rather than charge transfer, so the process 
will occur over a thicker region behind the shock. 
Thus O, N and especially Nc and He should ini- 
tially form ring distributions and be picked up by 
the plasma. As with the protons, the initial width 
of the ring varies as smO and the initial speed along 
the field as cosO. 

Heavy ions present in the upstream plasma can 
also have peculiar downstream distributions due 
to their passage through the shock. They are de- 
celerated by the electric potential jump associated 
with the shock, and because of their large mass to 
charge ratios they are decelerated less than the 



protons. Fusclier & Schmidt (1997) find that the 
initial ring velocity in this case is 

= VMm/q - 1) + l/16)/(m/g))i/2 (6) 

for a strong perpendicular shock. Thus we expect 
that heavy ions passing through the shock will 
have values of Vj_ between Vs and SVs /4. A few 
observations exist to test this expectation. The 
line widths of C IV and He II lines in the Hopkins 
Ultraviolet Telescope spectrum of SN1006 (Ray- 
mond et al. 1995) are the same as the width of 
Ha within substantial uncertainties, and the O VI 
line observed with FUSE is consistent with the 
same width (Korreck et al. 2004). Ghavamian et 
al. (2002) find that the pre-shock neutral fraction 
of H is about 0.1, while that of He is at least 0.7. 
Since C has a lower ionization potential than H, 
and O has the same ionization potential has H, 
these elements also have small pre-shock neutral 
fractions. Thus C and O should have larger val- 
ues of V±_ than H, while He should be primarily 
a pickup ion distribution. Unfortunately, the 10% 
to 30% uncertainties in the line widths preclude a 
definitive comparison, but with somewhat higher 
quality profiles for the UV lines one could begin 
to constrain the magnetic field direction. 

3.6. Cosmic Ray Modified Shocks 

Except for some consideration of magnetic field 

amplification, the discussion above assumes a sim- 
ple magnetohydrodynamic shock. However, both 
observations and theory (e.g., Drury & Yolk 1981; 
Malkov et al. 2000; Warren et al. 2005) indi- 
cate that a substantial fraction of the energy dis- 
sipated in the shock, as much as 80%, can be con- 
verted to cosmic rays. This results in a "modi- 
fied shock" structure with several interesting fea- 
tures; 1) a particle velocity distribution such as 
a Maxwellian with a power law tail, 2) a smooth 
transition rather than a sharp shock jump, 3) a 
compression ratio higher than the hydrodynamic 
factor of 4, and, 4) a lower proton temperature for 
a given shock speed, since less energy is available 
to heat the gas. 

The Balmcr line profiles arc not expected to re- 
veal the non-thermal tails predicted for strong dif- 
fusive shock acceleration of cosmic rays, since only 
a very small fraction of the particles (~ 10~^) are 
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accelerated. Also, the charge transfer cross sec- 
tion declines rapidly at speeds above about 2000 
km s^^ (e.g., Schultz et al. 2008), so the faster 
protons are less likely to produce broad compo- 
nent neutrals. Therefore, direct detection of the 
power law tail will be very difficult. 

The smooth transition could change the profile 
in a manner incompatible with observations, in 
that the gradually increasing temperature would 
give a composite Ha profile which is the sum of 
profiles formed at all the temperatures in the shock 
transition. It would probably not resemble the 
easily separable broad and narrow component pro- 
files observed. This difficulty would be avoided 
if the smooth transition occurs on a length scale 
smaller than the length scale for charge trans- 
fer, since few broad component neutrals would 
form in the intermediate temperature region. The 
length scale for a modified shock is k/Fs, where 
K is the cosmic ray diffusion coefficient. The 
charge transfer length scale is Vs/{npqcT), where 
qcT is the charge transfer rate coefficient. Since 
qcT ~ 3x 10"'' cm'^s"^ in the downstream plasma, 
K should be smaller than about 10^^ cm^s""'^. Val- 
ues of K of this order are required to accelerate 
cosmic rays to high energies within an SNR life- 
time, but they are comparable to the Bohm limit, 
and therefore at the low end of the range of plau- 
sible values. 

A high compression ratio, say 7 rather than the 
usual 4, would mean that the narrow component 
neutrals move at &Vs/1 relative to the postshock 
gas, rather than 3^5/4, so the PUI component 
will have larger initial parallel and perpendicular 
velocities by 14%. On the other hand, if a large 
fraction of the shock energy goes into nonthermal 
particles, the thermal speed of the protons will be 
reduced by a factor (1 -I- Pc/Pg)"^/^, where Pc 
and Pq are the cosmic ray and gas pressures. If 
Pc is comparable to Pq, the thermal part of the 
line width will be seriously affected and the shock 
speed will be underestimated if Pc is assumed to 
be small. Most of the Balmer line filaments stud- 
ied to date show very weak radio emission (e.g., 
the NW filament in SN1006 and the northern fila- 
ment in the Cygnus Loop; Ghavamian et al. 2001, 
2002), so Pc/Pg is probably small. 



4. Application to Tycho's SNR 

Tycho's supernova remnant presents a good 
opportunity to search for the effects described 
above becaus e of its relatively high pr e-shock neu- 
tral fraction ( Ghavamian et al. 2000l ). the excel- 
lent high and low resolution spectra of knot g 



(jSmith et al. 19911 : iGhavamian et al. 20011 : iLee et al. 20071 ). 

and the extensive X-ray and radio s tudies of both 

the thermal and non-thermal shocks (Pickcl fc Jones 19851: 



Dickel et al. 1991 



Hwang et al. 2002 : Warren et al. 2005 



Bamba et al. 2005[ ). The preshock density is ap- 
proximately 1 cm~'^, the pre-shock neutral frac- 
tion is approximately 0.85 and the shock speed is 
approximately 2000 km s~^ ( Ghavamian et al. 2000t 
[Gh avamian et al. 200ll) . The magnetic field is 
likely to be amplified in strong SNR shocks (e.g., 
Lucek & BeU 2000; Vink & Laming 2003), but its 
strength is not accurately known. Non-thermal 
synchrotron emission from nearby parts of the 
blast wave of Tycho's SNR, implies that the mag- 
netic field is on th e order of 100 jiG = IPioo 
(| Warren et al. 2005h . yielding an Alfven speed of 
approximately 100 km s~^. The field direction 
is not known with certainty, though Dickel et 
al. (1991) show that the polarization indicates 
a predominantly radial field on scales of a few arc- 
seconds behind the shock. The thermal pressure 
implied by the Rankine-Hugoniot jump conditions 
with the shock speed and the pre-shock density 
above yields a plasma /3 of 12/i?^QQ. 

There is a significant shift between the cen- 
troids of the broad and narrow components of the 
Ha profiles in Tycho. Smith et al. (1991) and 
Ghavamian et al. (2001) found shifts between the 
broad and narrow components of Ha of 240±60 
and 132±35 km s~^, respectively, for two slit po- 
sitions in knot g. Smith et al. interpreted the shift 
as an indication that the shock normal does not 
lie in the plane of the sky, so that the shift repre- 
sents a small component of the post-shock plasma 
speed. This interpretation is consistent with the 
observation of Lee et al. (2004), who showed that 
the centroid of the narrow component is shifted 
with respect to the centroid of the ambient gas 
in that region (though there is some uncertainty 
about the size of this shift; Lee et al. 2007). How- 
ever, the shock normal cannot be very far from the 
plane of the sky, since very strong limb brightening 
is required to account for the observed brightness 
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of knot g. 

Alternatively, it is possible that the shift be- 
tween broad and narrow component centroids is 
related to the projection of vtuik onto the line of 
sight. The shift is limited to approximately Va, 
so a shift of the magnitude measured would re- 
quire that the projection of the magnetic field di- 
rection onto the direction parallel the line of sight 
be fairly large, and therefore must be near 90°. 
Within the limits of the data now available, we 
cannot tell whether the shift between broad and 
narrow centroids is essentially a geometrical effect, 
as proposed by Smith et al. (1991) or a result of 
the pickup ion bulk speed discussed above. 

A second puzzle relates to the nature of turbu- 
lence downstream from the shock. If the 100 /iG 
field is generated by turbulent amplification in the 
shock front, it will be fairly disordered. The non- 
resonant mechanism proposed by Bell (2004) pre- 
dicts that the scale of the turbulence is smaller 
than the gyroradius of cosmic ray protons (Zi- 
rakashvili et al. 2008), and generally yields a 
perpendicular shock. Compression by the shock 
would also make the mean field direction more 
perpendicular to the shock normal. Giacalone & 
Jokipii (2007) and Zirakashvih & Ptuskin (2008) 
study the effects of density inhomogeneities on 
magnetic field generated downstream of the shock. 
Both works find significant magnetic amplifica- 
tion, and Zirakashvih & Ptuskin (2008) remark 
that the magnetic field component parallel to the 
shock normal is more enhanced. The interaction 
between the the pickup ions and the field also 
tends to bend the field toward the shock normal, 
and the observed field in Tycho's SNR is nearly 
radial at the edge of the SNR (Dickel et al. 1991.) . 
A turbulent field would suggest that the pickup 
process occurs over a large range of 0, smearing 
out the profile as in Figure 5d. Detection of a 
non-Gaussian profile in Ha would provide some 
idea of the nature of the turbulence. This will 
require very good data and careful assessment of 
the instrument profile and the background level, 
however, and existing data do not provide useful 
constraints. 



5. Discussion 

5.1. Caveats 

There are several qualifications to the analysis 
presented here. One is the use of PUI analysis 
based on Alfven waves, which is appropriate for 
a cold plasma. As mentioned above, (3 is around 
12 for Tycho's SNR, and that will be typical for 
the strong shocks seen as Balmer line filaments. 
Thus, other wave modes may be important. It is 
unknown whether they will tend to change the di- 
rections, rather than the energies, of the protons 
in the way that Alfven waves do. It is also pos- 
sible that they will provide better coupling to the 
electrons. 

Another question is whether the amplified 
B field behind the shock is strongly turbulent 
on small scales. If so, PUI w ould be gener- 
ated over a broad r ange of angles (jlsenberg 19981: 
Nemeth et al. 20"ool) , tending to wash out any line- 
shift signature in the Ha profile. The polarization 
measurements of Dickel et al (1991) indicate that 
the field is reasonably well ordered on the scale of 
their resolution, but it could be highly random on 
the 0.1" scale over which the Ha is produced. 

Finally, there is the question of momentum con- 
servation when a significant fraction of the down- 
stream plasma has been picked up and streams 
along an oblique magnetic field. In this case, the 
thermal plasma would presumably act to cancel 
the transverse momentum, resulting in a rotation 
of the field toward the shock normal. We plan to 
quantitatively investigate this interaction in the 
near future. The resolution may lie in the density 
gradient of the pickup ions, but further calcula- 
tions are needed. 

5.2. Implications for Balmer line filament 
analysis 

If the pickup ions provide a significant contri- 
bution to the Ha profile, values of Vs derived from 
Gaussian fits are somewhat in error. This error 
would propagate into distances derived by com- 
bining shock speeds derived from the Balmer line 
profiles with proper motions (e.g., Winkler, Gupta 
& Long 2003). The modifications are probably not 
severe, and in cases such as SN1006, where the pre- 
shock neutral fraction is only 10% and the contri- 
bution of pickup ions to the Balmer line profiles 



8 



is only 5%, they would be completely negligible. 
In cases where the pre-shock neutral fraction is of 
order 50%, as much as 25% of the broad compo- 
nent emission could arise from atoms produced by 
charge transfer from pickup ions. In such cases Vg 
would probably be underestimated. On the other 
hand, waves emitted by the pickup ions as they 
isotropize could heat the protons and lead to a 
compensating effect. 

If non-Maxwellian distributions can be ob- 
served by way of distortions of the Ha profiles of 
non-radiative shocks, they could contain unique 
information about the strength and direction of 
the magnetic field and the level of turbulence in 
the region where the Ha emission arises. The most 
promising SNR where non-Maxwellian distribu- 
tions might be found is probably Tycho, thanks 
to its large neutral fraction and relatively bright 
Ha emission. 
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Fig. 1. Bisphcrical distribution in velocity space 
where the shock speed equals 1.0. The Alfven 
speed for this example was 0.2 Vs and 9 was 60°. 

Fig. 2. — Bulk velocity of a bispherical distribu- 
tion along the field direction for a 2000 km s~^ 
shock with a pre-shock density of 1 cm^"' and a 
post-shock density of 4 cm^'' for a range of post- 
shock magnetic field strengths and angles between 
the shock normal and the field. The correspond- 
ing Alfven speeds are 880, 710, 530, 350, 180 and 
35 km s-^ 

Fig. 3. — Variation of proton and neutral hydro- 
gen densities behind a 2000 km s"'^ shock. This 
model does not include the effects of weighting 
with the charge transfer cross section or of drift on 
the pickup ions along the magnetic field, both of 
which tend to increase the velocity of the pickup 
ions relative to the shock front and reduce their 
density. 

Fig. 4. — Proton velocity distribution for an angle 
^=70° between the field and the shock normal, 
an Alfven speed Va=100 km s~^, and a pickup 
ion density of 0.25 the total density. The lower 
dashed line is the bispherical distribution, the up- 
per dashed line is the thermal proton distribution, 
and the solid line is the total. 

Fig. 5. — Velocity distributions for various com- 
binations of parameters, a) Ratios of pickup ions 
to thermal protons ranging from 0.05 (outermost 
curve)to 0.65, b) Angles between the field and the 
shock normal of 45° (outermost curve) and 75°, 
c) Alfven speeds of 0.05 Vs (outermost curve) and 
0.20 Vs, and d) the velocity distribution for a dis- 
tribution of angles assuming isotropic turbulence 
upstream and compression of B± by a factor of 4. 
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Fig. 6. — Energy of the bispherical distribution 

as a fraction of the initial energy of the neutral 
atoms. The energy is computed in the rest frame 
of the post-shock gas. The curves correspond to 
values of 9 of 85°, 75°, 65°, 55°, 45°, 35°, 25°, 15° 
and 5° from top to bottom. 
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